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The motivations for investigating a theory of gravitation based on a concept of "ether" 
are discussed- a crucial point is the existence of an alternative interpretation of special 
relativity, named the Lorentz-Poincare ether theory. The basic equations of one such 
theory of gravity, based on just one scalar field, are presented. To check this theory in 
celestial mechanics, an "asymptotic" scheme of post-Newtonian (PN) approximation is 
summarized and its difference with the standard PN scheme is emphasized. The deriva- 
tion of PN equations of motion for the mass centers, based on the asymptotic scheme, 
is outlined. They are implemented for the major bodies of the solar system and the 
prediction for Mercury is compared with an ephemeris based on general relativity. 
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1. Introduction: Why an Ether Theory of Gravitation? 

(i) Our first motivation was to extend the Lorentz-Poincare ether theory so that 
gravitation be included. The Lorentz-Poincare ether theory may be described as the 
theory according to which a) the ether is an inertial frame E such that Maxwell's 
equations are valid in E, and b) any material object that moves with respect to 
E undergoes a Lorentz contraction. As shown, in particular, by Prokhovnikjd this 
theory is physically equivalent to standard special relativity (SR). This makes its 
ether undetectable, hence "superfluous" (Einstein 1905). But since SR does not 
involve gravitation, we may ask whether a preferred frame could exist but remain 
hidden in the absence of gravitation, and become detectable in its presence. 

(ii) An "ether" could help to make quantum theory and gravitation theory com- 
patible. Quantum theory was originally built in a flat space-time, moreover it uses 
a preferred time. In flat space-time, this is the inertial time, which depends on the 
inertial frame in a way that remains compatible with Lorentz invariance. Yet in 
the curved space-time of gravitation, there seems to be no way to prefer some time 
coordinate, except if we a priori admit a "preferred space-time foliation" , i.e. an 
ether. This may be illustrated already for the case of the Klein-Gordon equation, 
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for which many possible extensions to curved space-time a priori exist, but one is 
preferred if we have an etherJa Further, quantum theory shows that "vacuum" has 
physical effects, e.g. the Casimir effect, now experimentally confirmed. 

(iii) Investigating a strongly alternative theory also opens a new way to solve 
some problems common to general relativity (GR) and to most extensions of it: 
a) Singularities: The investigated "scalar ether-theory" does avoid singularities, in 
gravitational collapseH and in cosmology as well.Q b) Gauge condition: The solutions 
to the underdeterminacy of the Einstein equations as a system of partial differen- 
tial equations are either to say that the Lorentz manifold is determined modulo 
diffeomorphisms, or to add a gauge condition in a fixed space-time manifold (the 
latter way is used in applications). What is the precise link between these two so- 
lutions? In the scalar ether-theory, there is no need for any gauge, yet space-time 
is fixed, c) Galactical dark matter: Identified candidates seem poorly found. In the 
present ether-theory, the preferred-frame effects are probably more important at 
the galactical scale and beyond, due to the large time scales involved. 

2. Basic Equations of the Investigated Theory 

N.B.: Most equations are preferred-frame ones with space covariance only. 

(i) Gravitation is seen as Archimedes' thrust in an imagined perfect fluid ("ether") 
with pressure p e and density p e — p e (p e ). This leads to define the gravity accelera- 
tion vector as follows :□ 

gradp e 

g = • (1) 

Pe 

Note that, due to Eq. (Q), p e and p e decrease towards the attraction, thus p e (x, t) < 
p^°(t) = Sup xGM/ 9 e (x, t) in a gravitational field, where M is the "space" manifold, 
i.e. the set of the positions x in the preferred frame E. 

(ii) Assumed metric effects of a gravitational /le/cf.i We assume that the space- 
time R x M is equipped with a flat metric 7 for which the preferred frame E is 
an inertial (Galilean) frame. The inertial time t in E is called the "absolute time" , 
and the Euclidean space metric associated with 7 in the frame E is denoted by 
g°. Yet we also assume that, in a gravitational field, i.e. p e (x,t) < pf(t), there 
are metric effects, similar to those due to uniform motion (see Ref. 1): the meters 
are contracted and the clocks are slowed down, in the ratio (3 = p e (x, t) / pf (t) and 
(for the meters) in the direction g only. This means a dilation (contraction) of the 
length (time) intervals, when they are indeed measured with physical instruments, 
as compared with those that would be evaluated in terms of the flat space metric 
g° and the absolute time t. Hence, the "physical" space metric g in the frame E is 
a Riemannian one, and the measured time is a "local" one, denoted by t x (at point 
x G M). Thus, the "physical" space-time metric 7 is a curved Lorentzian metric. 



6, 2008 4:57 WSPC/Guidclincs-IJMPA Joaotext 



Moreover, SR leads to assuming the relation p e = c p e . 

(iii) Gravitational field equation. The following equation is stated for the scalar 
gravitational field p e : 

A ^-^^=4^Gap e , (2) 

where a is the T 00 component of the energy-momentum tensor of matter and non- 
gravitational fields T, when the time coordinate is x° — ct with t the absolute time, 
and in any spatial coordinates that are adapted to the preferred reference frame 
E.B The derivative with respect to the local time is defined by 

dt x 0(t,x)df v ' 

In Eq. (|^), the Laplace (-Beltrami) operator is defined with the curved space metric 
g (relative to the frame E). The same is true for the grad operator in Eq. ([!]). 

(iv) Dynamics is governed by Newton's second law: force = time-derivative of 
momentum.! The force over the test particle is the gravitational force m(i>)g, plus 
the nongravitational (e.g. electromagnetic) force, where m(v) is the relativistic iner- 
tial mass, involving the Lorentz factor. The momentum is m{v )v; v and v = |v| are 
evaluated with the physical metric. The time-derivative of momentum is uniquelv 
defined from compelling requirements (including Leibniz' rule for a scalar product) .□ 
In the static case, that extension of Newton's 2nd law implies Einstein's geodesic 
motion. For a dust, we may apply this extension pointwise in the continuum, and 
it implies a new equation for continuum dynamics:! 

T; ; „ = b„, b ee \g jkfl Ti\ h ee - l -g lkfi T ok . (4) 

The universality of gravity is expressed in the fact that Eq. (^) is assumed to hold 
true for any material medium (thus also for a nongravitational field) . 

3. Asymptotic Post-Newtonian Approximation 

An "asymptotic" nost-Newtonian approximation (PNA) was developed! (cf. Fu- 
tamase & SchutzB and RendalH in GR; in GR, the local field equations of the 
asymptotic method have not been used to get equations of motion for the mass 
centers of extended bodies. This has beeen done in the present theoryEJ'EJ) : 

(i) The gravitational field and the matter fields are expanded (in the standard 
PNA,E3£j only the gravitational field is expanded). 



(ii) For definiteness, each body is assumed to be made of a barotropic perfect 
fluid (one fluid per body). Other constitutive laws may also be considered, of course. 
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(iii) A family (S A ) of gravitating systems is deduced from the given system 
S. To do this, we iise the fact that an exact similarity transformation exists in 
Newtonian gravity.Qtl This transformation is applied to the initial data for S. For 
this, V = c 2 (l — /)/2 is substituted for the Newtonian potential, where / = 700 
(1 - / < 1 for a weak field). The initial data for S (the system of interest, e.g. the 
solar system) is generaljj in contrast with Ref. 8 in which the initial space metric 
was very special. (In Ref. 9, the family was a priori assumed.) 

Adopting units [T] A = [Tj/A 1 / 2 and [M] A = A[M] for system S A , all fields are 
ord(A°), and the small parameter A is proportional to 1/c 2 (in fact A = (co/c) 2 , 



and [M]). Therefore, 
(In the standard PN 



where Cq is the velocity of light in the starting units [T . 
the derivation of asymptotic expansions is straightforward. t 
scheme, 1/c 2 is formally considered as a small parameter.) That 1/c 2 , not 1/c, turns 
out to be the effective small parameter, is due to the fact that it is only 1/c 2 that 
enters in the equations. The theory admits consistent expansions in powers of A (or 
1/c 2 ). The first (zero-order) term is Newtonian gravity, hence the theory admits a 
correct Newtonian limit. The first-order approximation in A or 1/c 2 is the first PN 
approximation. Using these expansions is justified insofar as the system of interest 
corresponds to a small value Aq of A, which is the case, e.g., for the solar system. 



4. PN Equations of Motion for the Mass Centers 

The mass centers (MC) are defined^ as local barycenters of the rest-mass density 
Pexact (instead of, e.g., the active energy density a, Eq. (§)), because (i) p e xact = 
outside the bodies (which is wrong if one takes instead a density that involves 
gravitational energy, since the latter is distributed in the whole space) and (ii) Pexact, 
or rather its PN approximation, obeys the usual continuity equation, i.e. without 
adding gravitational energy and its flux. To get the MC's equation of motion, one 
just integrates the local equations of motion inside the different bodiesJlIM Due 
to the use of the asymptotic method, the local equations of orders and 1 in A or 
1/c 2 are separated. E.g.: 

d t po + dj{p u 3 ) = 0, d t pi + dj{pxu{ + p Q u{) = (5) 

for the continuity equation, derived from the time component at the first PNA. The 
equation of order 1 is linear with respect to the fields of order 1. Therefore, separate 
equations are also obtained for the MC's, and the equation for PN corrections (order 
1) is linear with respect to order-1 quantities. To get tractable equations, every field 
is decomposed into a self-field and an external field, and account is taken of the 
"good separation" between different bodies, which means that 

n = Sup a ^ b (r a /|x a - x b |) < 1 (6) 

(r a is the radius of body (a)). In the solar system, terms up to and including rf 
must be retained.Eil A rigid motion, possibly including self-rotation, is assumed for 
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each body. Finally, the rest-mass density of the order 0, po, is assumed spherical 
for each body, at the stage of calculating the PN corrections. One thus gets explicit 
equations of motion for the mass centers.!!!! They show that the self-rotation of the 
bodies and their internal structure influence the motion from the first PNA. This 
follows naturally from using the asymptotic method and should hold true for GR. 

5. Implementation. Comparison with a General-Relativistic 
Ephemeris 

In order to use the equations of motion for the mass centers0 so as to check the 
theory, we have to know the values of the parameters that enter these equations. 
These are the 0-order masses M a of the bodies (here the major bodies of the solar 
system), the initial conditions of their motion, and the constant velocity V of the 
global zero-order mass center of the solar system, with respect to the preferred frame 
E (and also the constant G).0 (Of course there is no parameter like V in conven- 
tional theories.) These unknown parameters depend on the theory. They must be 
determined by optimizing the agreement between predictions and observations.!!!!! 
Our computer code loops on the numerical solution of the translational equations 
of motion in order to optimize the parameters-El This code has been tested by in- 
vestigating in which measure one may reproduce (over one century) the predictions 
of the DE403 ephemerisjil by using purely Newtonian equations of motionJll It has 
also been applied to adjust over 60 centuries a less simplified model, in jazhich the 
PN corrections in the Schwarzschild field of the Sun are also considered.t3 

In the version of the code that incorporates the equations of motionE derived 
from the present theory, a Lorentz transform allows to pass from the preferred 
reference frame to the frame bound with the zero-order global barycenter, and 
vice-versa. This transform, as well as the inverse transform, is determined by the 
adjustable vector V. Thus, the adjustment process of the translational equations 
on observational data provides us eventually with the value of V that minimizes 
the residual with the set of observations. Note that the "observational data" are 
currently taken from an ephemeris based on GRp3 specifically we take a set of 
heliocentric positions of the eight major planets, between 1956 and 2000. With these 
input data, themselves a fitting of observations by GR equations, the magnitude of 
the optimal vector V is |V| ~ 3 km/s, which is significant. The difference between 
DE403 and our thus-adjusted equations of motion is small (cf. the residual advance 
in Mercury's longitude of perihelion, with respect to Newton's theory: 43"), but 
significant (Fig. 1. The self-rotation of all nine bodies is neglected). Our current 
project is to adjust the theory on a set of true astronomical observations. 
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Fig. 1. Difference between the scalar theory and the DE403 cphemeris of the JPL. 
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